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lO ■ Abstract 



We have surveyed spatial profiles of the Fe Ka lines in the Galactic center diffuse X-rays (GCDX), 

including the transient region from the GCDX to the Galactic ridge X-ray emission (GRXE), with the 

[^ ' Suzaku satellite. We resolved Fc Ka line complex into three lines of Fel, FeXXV and FeXXVI Ka, and 

obtained their spatial intensity profiles with the resolution of ~ 0.°1. We compared the FeXXV Ka profile 

with a stellar mass distribution (SMD) model made from near infrared observations. The intensity profile 

^ '. of FeXXV Ka is nicely fitted with the SMD model in the GRXE region, while that in the GCDX region 

I ■ shows 3.8±0.3 (0.°2 < |;| < l.°5) or 19±6 {\l\ < 0.°2) times excess over the best-fit SMD model in the GRXE 

2 I region. Thus FeXXV Ka in the GCDX is hardly explained by the same origin of the GRXE. In the case of 

*-^ . point source origin, a new population with the extremely strong FeXXV Ka line is required. An alternative 

jrt ' possibility is that the majority of the GCDX is truly diffuse optically thin thermal plasma. 

Key words: Galaxy: Center — Inter stellar medium — X-ray spectra 



temperature variations and/or may be contaminated by 
Diffuse X-rays were discovered from the inner disk of the possible non-thermal emissions, we concentrate on the flux 
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<*' ' 1. Introduction Since spectra of the GCDX and GRXE have spatial 

oo ■ 

o 

00 

p.^ ,. Milky Way Galaxy with the HEAO-1 satellite (Worrall ct of FeXXV Ka line as a pure component of the thin ther- 

al. 1982), and called as the Galactic ridge X-ray emis- mal plasma. Previous studies of the Fe Ka lines with 

sion (GRXE). EXOSAT observed more details and found the Ginga and RXTE satellites (Yamauchi & Koyama 

that the GRXE were extended by about ±60° in longi- 1993; Revnivtsev et al. 2006) had limited spatial reso- 

tude and ±(l°-2°) in latitude (Warwick et al. 1985). The lutions (FWHM ~ l°-2°) compared to the size of the 

Tenma satellite discovered strong emission lines at the en- GCDX, and the FeXXV Ka line flux was contaminated 

ergy of about 6.7 kcV (here the 6. 7-kcV line) in the GRXE by the adjacent Fel and XXVI Ka lines due to the lim- 

S^ i (Koyama et al. 1986). The Ginga satellite made survey ited energy resolutions (FWHM ~ 1.1 keV at 6 kcV). 

;_i ' observations along the Galactic plane with the 6.7-keV Yamauchi et al. (2009) reported that the intensity ratio of 

5t 1 line, and found a flux peak at the Galactic center with the the three Fe Ka lines has spatial dependence and hence 

width of about 1°. We call this peak as the Galactic center to resolve the three lines is essential for the study the 

diffuse X-rays (GCDX) (Koyama et al. 1989; Yamauchi et nature of FeXXV Ka line. This paper reports, for the 

al. 1990; Yamauchi & Koyama 1993). first time, the FeXXV Ka line distributions in the GCDX 

The ASCA satellite separated the 6.7-kcV line in the and the transient region from the GCDX to the GRXE 

GCDX into three emission fines at 6.4, 6.7 and 7.0 keV, (-3° <l <2°, -2° < 6 < 1°) with the CCD cameras on- 

which were respectively identified as Ka lines from neu- board the Suzaku satellite. 

tral (Fel), He-like (FeXXV) and H-like (FeXXVl) iron The observations and data reduction are described in 
(Koyama et al. 1996). For the GRXE, the Suzaku satel- section 2. The analysis and the intensities of Fel Ka, 
lite resolved the Fe Ka lines (Ebisawa et al 2008). The FeXXV Ka and FeXXVI Ka are given in section 3. In sec- 
strongest emission line in the GCDX and GRXE is FeXXV tion 4, we compare the FeXXV Ka profiles of the GCDX 
Ka at 6.7 kcV, and hence would be due to an optically and GRXE with that of a stellar mass model made from 
thin thermal plasma. From the fiux ratio of FeXXVI near infrared (NIR) observations. Based on the profile 
Ka/FcXXV Ka, the plasma temperature of the GCDX comparison between NIR stars and FeXXV Ka, we dis- 
was estimated to be around kT =5-10 keV (e.g. Koyama cuss origin of the GCDX. In this paper, uncertainties 
et al. 2007b; Yamauchi et al. 2009). are quoted at the 90% confidence range unless otherwise 
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stated. 

2. Observations and Data Reduction 

The data in this paper were obtained from September 

2005 to March 2009 with the X-ray Imaging Spectrometer 
(XIS) onboard the Suzaku satellite. Most of the data were 
observed as the Suzaku Key project observation. The to- 
tal effective exposure is about 3.2 Ms. The details of these 
observations are summarized in table 1. 

The XIS contains four sets of X-ray CCD camera sys- 
tems (XIS 0, 1, 2, and 3) placed on the focal planes of 
four X-Ray Telescopes (XRTs) onboard the Suzaku satel- 
lite. XIS 0, 2, and 3 have front-illuminated (FI) CCDs, 
while XIS 1 has a back-illuminated (BI) CCD. Detailed 
descriptions of the Suzaku satellite, the XRT, and the 
XIS can be found in Mitsuda et al. (2007), Serlemitsos et 
al. (2007), and Koyama et al. (2007a), respectively. 

In these observations, the XIS was operated in the nor- 
mal clocking full window mode, with the time resolution 
of 8.0 s. One of the FI CCD cameras (XIS 2) has been out 
of function since November 2006, and hence 71% of the 
total exposure was observed without XIS 2. Although the 
XIS CCDs were significantly degraded by on-orbit parti- 
cle radiation, the CCD performance had been calibrated 
well with the checker-flag charge-injection (CI) method 
(Nakajima et al. 2008; Ozawa et al. 2009) until October 

2006 , and has been restored with the spaced-row CI (SCI) 
method (Prigozhin et al. 2008; Uchiyama et al. 2009). The 
SCI method was used in 75% of the total exposure. The 
XIS 2 and SCI status in the observations is shown in table 
1. 

We reprocessed the data using the calibration data 
base (CALDB) released at 2008-10-20 and xispi in the 
HE AD AS software package version 6.6.1. We confirmed 
that the center energy of Mn I Ka line of the onboard cali- 
bration sources (^^Fe) was consistent with 5.895 eV within 
the calibration uncertainty of 20 eV, and the energy reso- 
lutions were 130-180 eV (FWHM) at 5.9 keV during the 
observations. 

3. Analysis and Results 

In order to investigate intensity profiles of Fe Ka lines 
along the Galactic longitude and latitude. X-ray spec- 
tra were made from the rectangles shown in figures la 
and b, respectively. The sizes of the rectangles along the 
longitude (figure la) are 6' x 12' (0.°1 in longitude and 
0.°2 in latitude), while those along the latitude (figure 
lb) are 12' x 6' (0.°2 in longitude and 0.°1 in latitude). 
Here we adopted a new coordinate as (/*,6») = (/ + 0.°056, 
&-|-0.°046), where the origin of this coordinate (/,,fe,) = 
(0°,0°) is the position of Sgr A* (Reid & Brunthaler 2004) 
and the 6* = 0° line corresponds the Galactic plane. With 
this new coordinate, the centers of each rectangle along 
the longitude (figure la) are given as (Z*,6*) = (0.°1 x n, 
0°) where n is an integer number. The centers of each 
rectangle along the latitude (figure lb) are &, = 0.°1 x n 
and U = -l.°104, -0.°784, -0.°504, -0.°114, 0.°226, 0.°556, 



0.°876 and l.°226. To extend the profile to larger 6* values, 
we used the data at ih,K.) = (0.°056, -l.°954) (OBSID 
502059010). Since the Fe Ka fiux at this position is lower 
than any other regions, we made a spectrum of 17' x 17' 
rectangle. To make spectra of diffuse emission, bright 
sources listed in table 2 and the CCD corners illuminated 
by the calibration sources were excluded from all the rect- 
angles. 
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Fig. 1. (a) The rectangles to extract X-ray spectra for the 
longitudinal distribution overlaid on the Suzaku XIS image in 
the 6.55-6.8 keV band. Bright point source regions (circles, 
see table 2) are excluded before making the spectra, (b) Same 
as (a) but for the latitudinal distribution. 

Using xissimarf gen (Ishisaki et al. 2007) and 
xisrmf gen, the effective area of the XRT and the response 
of the XIS were calculated, and ancillary response files and 
redistribution matrix files were made for each observation. 
Since the responses of the FIs are almost the same, we 
merged the FI spectra. 

The non X-ray background (NXB) spectra sorted by 
the cut-off rigidity were made from the same source region 
in the detector coordinate using xisnxbgen (Tawa et al. 
2008). The spectra of NXB were subtracted from those of 
the GCDX. 

We fitted the spectra in the 5-10 keV band with a phe- 
nomenological model which consists of four Gaussians and 
an absorbed power-law continuum. The line centers of 
the Gaussians were fixed to 6.40 (Fel Ka), 6.68 (FeXXV 
Ka), 6.97 (FeXXVI Ka) and 7.06 (Fel K/3) keV according 
to the APEC model (Smith et al. 2001) and Kaastra & 
Mewe (1993). The intrinsic widths of these lines were fixed 
to eV. The intensities of the Fel Ka, FeXXV Ka, and 
FeXXVI Ka lines in the model were free parameters. The 
intensity of the Fel K/3 line was fixed to 0.125 times that 
of the Fel Ka line according to Kaastra & Mewe (1993). 
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Table 1. Observation data list. 



OBSID 




Pointing direction 




Observation 


E.xp. 


XIS2* 


scr 


Object name 




<^2oao.o{°) 


*2000.o(°) 


in 


bn 


Start (UT) 


End (UT) 


',¥', 








1000:^7010 


266.51 


-iHM 


0.05 


-0.07 


2005-0<J-23T07:18:2.5 


»05-09-MTll:05:19 


44.7 


ok 


oil' 


ygr A East 


100027020 


266.30 


-29.17 


-0.25 


-0.05 


2005-09-24T14:17:17 


2005-09-25T17:27:19 


37.4 


ok 


off 


Sgr A west 


100037010 


266.30 


-29.17 


-0.25 


-0.05 


2005-09-29T04:35:41 


2005-09-30T04:29:19 


43.5 


ok 


off 


Sgr A west 


100037030 


266.20 


-29.35 


-0.45 


-0.07 


2005-09-30T06:06:32 


2005-09-30T07:41:20 


3.0 


ok 


off 


KS1741-293(GCBGD2) 


100037040 


266.51 


-28.93 


0.06 


-0.07 


2005-09-30T07:43:01 


2005-10-01T06:21:24 


42.9 


ok 


off 


Sgr A East 


100037060 


266.88 


-28.44 


0.64 


-0.10 


2005-10-10T12:28:01 


2005-10-12T07:05:23 


76.4 


ok 


off 


Sgr B2 


500018010 


266.16 


-29.47 


-0.57 


-0.09 


2006-02-20T12:45:25 


2006-02-23T10:50:14 


106.9 


ok 


off 


SGRC 


500005010 


266.77 


-28.63 


0.43 


-0.12 


2006-03-27T23:00:22 


2006-03-29T18:12:15 


88.4 


ok 


off 


GC CENTER #2 


100048010 


266.51 


-28.93 


0.06 


-0.07 


2006-09-08T02:23:24 


2006-09-09T09:06:15 


63.0 


ok 


off 


GALACTIC CENTER 


501040010 


266.69 


-28.38 


0.61 


0.07 


2006-09-21T17:29:01 


2006-09-23T06:52:20 


61.4 


ok 


on 


GC SGR B NORTH 


501040020 


266.69 


-28.38 


0.61 


0.07 


2006-09-24T05:03:12 


2006-09-25T07:00:24 


44.8 


ok 


on 


GC SGR B NORTH 


501008010 


266.50 


-29.17 


-0.15 


-0.19 


2006-09-26T14:18:16 


2006-09-29T21:25:14 


129.6 


ok 


on 


GC SOUTH 


501009010 


266.19 


-28.91 


-0.07 


0.18 


2006-09-29T21:26:07 


2006-10-01T06:55:19 


51.2 


ok 


on 


GC SOUTH BCD 


501010010 


266.26 


-30.37 


-1.29 


-0.64 


2006-10-07T02:16:52 


2006-10-08X10:19:19 


50.7 


ok 


off 


HESS J1745-303 


501049010 


265.38 


-29.75 


-1.17 


0.33 


2006-10-08T10:22:40 


2006-10-09T02:19:24 


19.6 


ok 


on 


GALACTIC CENTER 


501050010 


265.91 


-29.65 


-0.83 


0.00 


2006-10-09T02:20:25 


2006-10-09T13:39:24 


22.0 


ok 


on 


GALACTIC CENTER 


501051010 


265.70 


-29.93 


-1.17 


-0.00 


2006-10-09T13:40:09 


2006-10-10T06:44:24 


21.9 


ok 


on 


GALACTIC CENTER 


501052010 


265.50 


-30.21 


-1.50 


0.00 


2006-10-10T06:45:09 


2006-10-10T21:18:14 


19.3 


ok 


on 


GALACTIC CENTER 


501053010 


265.30 


-30..50 


-1.83 


-0.00 


2006-10-10T21:18:59 


2006-10-llT10:06:14 


21.9 


ok 


on 


GALACTIC CENTER 


501057010 


266.03 


-30.11 


-1.17 


-0.33 


2006-10-llT10:07:27 


2006-10-12T03:28:14 


20.5 


ok 


on 


GALACTIC CENTER 


501039010 


267.02 


-28.35 


0.78 


-0.16 


2007-03-03T12:20:20 


2007-03-05T12:39:25 


96.4 




on 


GC SGR B EAST 


501046010 


265.98 


-28.90 


-0.17 


0.33 


2007-03-10T15:03:10 


2007-03-llT03:55:14 


25.2 




on 


GALACTIC CENTER 


501047010 


265.78 


-29.19 


-0.50 


0.33 


2007-03-llT03:55:59 


2007-03-llT19:04:14 


25.6 




on 


GALACTIC CENTER GC2 


501048010 


265.58 


-29.47 


-0.83 


0.33 


2007-03-llT19:04:59 


2007-03-12T08:09:14 


27.5 




on 


GALACTIC CENTER GC3 


501054010 


266.63 


-29.25 


-0.17 


-0.33 


2007-03-12T08:ll:07 


2007-03-12T23:58:24 


26.1 




on 


GALACTIC CENTER GC9 


501055010 


266.43 


-29.-54 


-0.50 


-0.33 


2007-03-12T23:59:09 


2007-03-13T15:40:19 


27.2 




on 


GALACTIC CENTER GCIO 


501056010 


266.23 


-29.82 


-0.83 


-0.33 


2007-03-13T15:41:12 


2007-03-14T05:00:24 


26.5 




on 


GALACTIC CENTER GCll 


501058010 


266.98 


-27.72 


1.30 


0.20 


2007-03-14T05:02:29 


2007-03-15X18:55:14 


63.3 




on 


GC SGR D NORTH 


501059010 


267.09 


-27.94 


1.17 


0.00 


2007-03-15T18:55:51 


2007-03-17T05:06:19 


62.2 




on 


GC SGR D 


501060010 


267.29 


-27.65 


1.50 


0.00 


2007-03-17T05:07:04 


2007-03-18T20:58:14 


64.8 




on 


GC SGR D EAST 


502022010 


266.81 


-28.88 


0.23 


-0.27 


2007-08-31T12:33:33 


2007-09-03T19:00:25 


134.8 




on 


{L,B)={0. 25,-0. 27) 


102013010 


266.51 


-28.93 


0.06 


-0.07 


2007-09-03T19:01:10 


2007-09-05T05:20:20 


51.4 




on 


GALACTIC CENTER 


502020010 


267.19 


-28.13 


1.05 


-0.17 


2007-09-06TOO:26:47 


2007-09-09T16:10:19 


139.1 




on 


SGR D SNR 


502002010 


267.16 


-29.14 


0.17 


-0.67 


2007-10-09T16:40:54 


2007-10-10T03:40:24 


23.2 




on 


GC14 


502003010 


266.96 


-29.42 


-0.17 


-0.67 


2007-10-10T03:41;13 


2007-10-10T15:20:24 


21.5 




on 


GC15 


502004010 


267.48 


-29.31 


0.17 


-1.00 


2007-10-10T15:21:17 


2007-10-llT01:00:24 


19.9 




on 


GC16 


502005010 


267.29 


-29.60 


-0.17 


-1.00 


2007-10-llT01:01:17 


2007-10-llTll:32:20 


20.6 




on 


GC17 


502006010 


266.18 


-28.62 


0.17 


0.33 


2007-10-llTll:34:01 


2007-10-llT23:07:14 


22.6 




on 


GC18 


502007010 


265.86 


-28.45 


0.17 


0.67 


2007-10-llT23:09:15 


2007-10-12T09:52:14 


22.0 




on 


GC19 


502008010 


265.66 


-28.73 


-0.17 


0.67 


2007-10-12T09:52;59 


2007-10-12T21:50:19 


23.8 




on 


GC20 


502009010 


267.48 


-27.37 


1.83 


-0.00 


2007-10-12T21:52:24 


2007-10-13T07:30:19 


20.9 




on 


GC21 


502010010 


266.38 


-28.34 


0.50 


0.33 


2007-10-13T07:32:00 


2007-10-13T18:50:24 


21.6 




on 


GC22 


502011010 


266.57 


-28.05 


0.83 


0.33 


2007-10-13T18:51:09 


2007-10-14T05:30:24 


23.0 




on 


GC23 


502016010 


266.23 


-30.11 


-1.08 


-0.48 


2008-03-02T18:08:00 


2008-03-04X17:40:19 


70.5 




on 


HESS J1745-303 1 


502017010 


266.47 


-30.09 


-0.95 


-0.65 


2008-03-06T13:26:36 


2008-03-08T16:00:24 


72.6 




on 


HESS J1745-303 2 


502018010 


266.06 


-30.24 


-1.27 


-0.43 


2008-03-08T16:02:17 


2008-03-10T21:00:19 


79.0 




on 


HESS J1745-303 3 


502051010 


266.94 


-28.15 


0.91 


0.01 


2008-03-11T06:19:45 


2008-03-15T05:30:18 


138.8 




on 


GC G0.9-t0.1 


502059010 


268.38 


-29.96 


-0.00 


-2.00 


2007-09-29T01:40:51 


2007-10-02X14:10:16 


136.8 




on 


GALACTIC BULGE 


503007010 


266.44 


-28.57 


0.33 


0.17 


2008-09-02T10:15:27 


2008-09-03T22:52:24 


52.2 




on 


GC LARGEPROJECTl 


503008010 


266.78 


-29.13 


0.00 


-0.38 


2008-09-03T22:53:29 


2008-09-05T06:56:19 


53.7 




on 


GC LARGEPR0.IECT2 


503009010 


266.45 


-29.34 


-0.33 


-0.24 


2008-09-05T06:57:08 


2008-09-06T15:55:24 


52.4 




on 


GC LARGEPR0.IECT3 


503010010 


266.04 


-29.55 


-0.70 


-0.05 


2008-09-06T15:56;13 


2008-09-08T01:39:24 


53.1 




on 


GC LARGEPR0.1ECT4 


503011010 


265.95 


-29.83 


-0.97 


-0.13 


2008-09-08T09:08:09 


2008-09-09T21:33:19 


57.6 




on 


GC LARGEPR0.IECT5 


503012010 


266.30 


-29.94 


-0.91 


-0.44 


2008-09-14T19:35:07 


2008-09-16T00:50:14 


57.7 




on 


GC LARGEPR0.IECT6 


503013010 


265.67 


-30.07 


-1.30 


-0.05 


2008-09-16X00:51:19 


2008-09-18T04:44:24 


104.8 




on 


GC LARGEPR0.IECT7 


503014010 


265.18 


-30.75 


-2.10 


-0.05 


2008-09-18T04:46:49 


2008-09-19T07:32:20 


55.4 




on 


GC LARGEPR0.IECT8 


503015010 


265.03 


-30.96 


-2.35 


-0.05 


2008-09-19X07:33:05 


2008-09-20T09:56:13 


56.8 




on 


GC LARGEPR0.IECT9 


503016010 


264.87 


-31.17 


-2.60 


-0.05 


2008-09-22T06:47:49 


2008-09-23T08:07:17 


52.2 




on 


GC LARGEPROJECTIO 


503017010 


264.72 


-31.38 


-2.85 


-0.05 


2008-09-23T08:08:10 


2008-09-24T09:21:13 


51.3 




on 


GC LARGEPROJECTl 1 


503018010 


264.56 


-31.60 


-3.10 


-0.05 


2008-09-24T09:27:54 


2008-09-24T22:30:24 


29.4 




on 


GC LARGEPR0JECT12 


503072010 


265.99 


-29.21 


-0.42 


0.17 


2009-03-06T02:39:11 


2009-03-09T02:55:24 


140.6 




on 


EXTENDED CHIMNEY 


403009010 


266.48 


-28.78 


0.17 


0.02 


2009-03-21T02:03:27 


2009-03-23T18:26:09 


110.8 




on 


ARCHES CLUSTER 



'^The XIS2 column sliows whether the XIS 2 worked 
^Tlie SCI column sliows whether the SCI mode was 



( "ok" ) or not in the observation, 
on or off in the observation. 



The normalization, photon index, and interstellar absorp- 
tion of the absorbed power-law model were free. The cross 
section of the photoelectric absorption was obtained from 
Morrison & McCammon (1983). 

The FIs and BI spectra were fitted simultaneously. We 
did not use the BI data above 8 keV, because the NXB 
fluxes become high in this energy band. If one rectangle 
region overlaps more than two separate observations, we 
made separate spectra and were fitted simultaneously. 

Some examples of the spectra and the best-fit models 
are shown in figure 2. In these figures, we find some resid- 
uals at 8-9 keV band, which are identified as FeXXV K/3, 
Ni XXVII Ka and other weak lines (Koyama et al. 2007b). 
Inclusion of these lines, however, does not change the best- 
fit intensities of the Fe Ka lines. The best-fit intensities of 
the three iron lines and the 5-10 keV band are shown in 
table 3. The profiles of the Fe XXV Ka line intensity along 



the Galactic longitude (6* = 0°) and latitude (/* = -0.°114) 
are shown in figure 3 as examples. 

We see global symmetry between the intensities in the 
positive and negative U (figure 3a) and 6* (figure 3b), 
except for the region of |/*| = 0.°l-0.°3 in figure 3a (see 
section 4.3). 

4. Discussion 

4-1. The distributions of the Fe XXV Ka line and the stel- 
lar mass 

Using the intensity profile, origin of the FeXXV Ka line 
in the GCDX is discussed in comparison with that in the 
GRXE. We fitted the profile of the FeXXV Ka hnc with 
a stellar mass distribution (SMD) model which was orig- 
inally compiled by Muno et al. (2006), based on the NIR 
observations (Launhardt et al. 2002; Kent et al. 1991). 
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Table 2. The positions and radii of bright X-ray sources excluded before making the diffuse X-ray spectra. 
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Name 


Position 
I b 


Radius 


IE 1740.7-2942 


-0.°887 


-0.°097 


4' 


CXOGCS J174445. 5-295042 


-0.°774 


-0.°062 


2' 


2E 1742.9-2929 


-0.°439 


-0.°383 


5' 


KS 1741-293 


-0.°432 


-0.°074 


3' 


Sgr A East 


-0.°047 


-0.°053 


3' 


Arches cluster 


0.°123 


0.°024 


2' 


IE 1743.1-2843 


0.°256 


-0.°029 


3' 


G0.61-f-0.01 


0.°620 


0.°016 


4' 


CXOGC J174645.3-281546 


0.°712 


0.°141 


2' 


HD 161507 


0.°837 


0.°085 


2' 


G0.9-h0.1 


0.°866 


0.°077 


2' 


SAX J1748.2-2808 


0.°994 


-0.°086 


2' 


AX J1749. 1-2733 


l.°592 


-0.°062 


2' 



tt,W=(l°-6OO,0°.O00) 



(Z„fe.)=(l°-226,-0°.200) 




S 0.01 r 




Energy (keV) 



Energy (keV) 



Fig. 2. Left panel: an example of the spectra for the longitudinal distribution with the best-fit model (solid line). Dotted Gaussians 
are Ka lines of Fel, Fexxv and Fexxvi. The weak line in the right is Pel K/3. For simplicity, only the FIs spectra are shown. The 
vertical bars show the Icr errors. Residuals from the best-fit model near at 8-9 keV are due to Nixxvil Ka and Fexxv K/3 and 
other weak lines (see text). Right panel: the same as the left panel, but that for the latitudinal distribution data. 

Table 3. List of observed intensities.* 



Galactic coordinate 



FelKa 



Observed intensity 
Fe XXV Ka 



^ 

Fe XXVI Ka 5-10 kcV 



Flag^ 



-0.156 


-0.046 


-0.100 


0.000 


7.31 


(6.83-7.72) 


18.4 


(17.7-18.9) 


7.97 


(7.44-8.42) 


20.7 


I 


-0.056 


-0.046 


0.000 


0.000 


12.4 


(11.5-13.3) 


45.5 


(44.3-46.7) 


13.4 


(12.9-14.3) 


44.3 


I 


0.044 


-0.046 


0.100 


0.000 


35.2 


(34.8-35.7) 


24.3 


(23.9-24.7) 


9.49 


(9.14-9.85) 


39.0 


I 


0.500 


-0.146 


0.556 


-0.100 


8.35 


(8.05-8.65) 


8.05 


(7.75-8.36) 


3.14 


(2.88-3.41) 


11.3 


b 


0.500 


-0.046 


0.556 


0.000 


9.95 


(9.71-10.3) 


10.2 


(9.95-10.5) 


3.25 


(3.01-3.52) 


13.7 


b 


0.500 


0.054 


0.556 


0.100 


4.79 


(4.51-5.02) 


9.06 


(8.68-9.32) 


2.28 


(2.02-2.53) 


9.83 


b 


0.000 


-2.000 


0.056 


-1.954 


0.04 


(0.00-0.12) 


0.64 


(0.55-0.74) 


0.27 


(0.18-0.36) 


1.08 


o 



*This table is available in its entirety as a machine-readable form in the online journal (see the PAS J web site). A 

portion is shown here for guidance regarding its form and contents. 

'I'The interstellar absorption is not corrected. The units are 10~^ photons s~^ cm~^ arcmin^^ and 10^^"* erg s^^ cm^^ 

arcmin^^, for the Fe Ka lines and the 5-10 keV band, respectively. Parentheses show the 90% error ranges. 

■I- The flags I , b and o mean that the data arc obtained along the Galactic longitude (figure la), the latitude (figure 

lb) and from {l^X) = (0.°056, -l.°954), respectively. 
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«••«» 




0.1 1 

Distance from Sgr A along the longitude U (degree) 
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0.1 1 
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Fig. 3. (a) Observed intensity profiles of the Fcxxv K« along the Galactic longitude (b» = 0°). The bars show the Icr errors, (b) 
The same as (a) but for the latitude distribution at Z* = — 0.°114. 



The SMD model consists of four components of the nu- 
clear stellar cluster (NSC), the nuclear stellar disk (NSD), 
the Galactic bulge (GB) and the Galactic disk (CD). The 
details of the SMD model arc described in Appendix. 

Wc made simultaneous fit for the longitudinal distribu- 
tion (figure la) and for many latitudinal distributions (fig- 
ure lb). In this fitting, we included the results of Koyama 
etal. (2007b) (|;*| <0.°4), Yamauchiet al. (2009) (1^*1 > 8°) 
and the result of {h,K) = (0.°056, -l.°954). These data 
were added in the longitudinal plot in figures 4a and c, 
except for those of Yamauchi ct al. (2009) of |6| > 0.°15. 
The additional data points in figures 4a and c are not ex- 
actly on 6* = 0°. However, in the simultaneous fitting, the 
real coordinate values were taken into account. 
4.1.1. The GRXE region 

The intensity distribution of the Fe Ka line along the 
Galactic plane has been already reported by Yamauchi 
& Koyama (1993) and Rcvnivtsev et al. (2006). Since 
Revnivtsev et al. (2006) found that the Fe Ka line inten- 
sity profile agrees well with the NIR profile in the GRXE 
region using the RXTE data, we first tried to fit the data 
of |Z*| > l.°5 regions and |6*| > 0.°5 regions simultaneously 
with the SMD model. We assumed that the mass emissiv- 
ities (see equation A7) for each component are the same, 
i.e. £nsc = £nsd = £gd = £gb = £■ Thus the free param- 
eter was only one, e or the normalization. The best-fit 
parameter was e = (5.85 ± 0.03) x lO'^'^ photons s~^ Mq^ . 
The fitting result for the longitudinal distribution is given 
in figure 4a, while that of the latitude at U = — 0.°114 is 
given in figure 4b as a typical example. 

As is shown in figure 4a, we obtained a nice fit in the 
GRXE data of |;*| > l.°5, in which the CD and partly GB 
components mainly contribute in the SMD model. Figure 
4b also shows reasonable fit in the region of |6*| > 0.°6, 
where the CD and GB components contribute equally. 

We checked the consistency between the Suzaku results 
(figure 4) and those of Revnivtsev et al. (2006) (figure 3 in 
their paper). The data in the GRXE region of |^*| = 10°- 
20° were compared since the limited spatial resolution of 
RXTE made it impossible to resolve the detailed structure 



around the Galactic center (GCDX). The result of figure 
4a is purely for FeXXV Ka line and the FeXXV Ka line 
intensity at |Z*| = 10°-20° is about 1 x 10~^ photons cm~^ 
s~^ arcmin"^, while that of Revnivtsev et al. (2006) is the 
sum of FcI,XXV and XXVI Ka lines of about 4.5 x 10"'' 
photons cm~^ s~^ dcg ~^ (hence 1.3 x 10"^ photons cm"^ 
s~^ arcmin^^). The Fe Ka line complex in the RXTE 
data was not separated into the three lines because of the 
limited energy resolution. From table 1 of Ebisawa et al 
(2008) and also from figure 3 of Yamauchi et al. (2009), 
we estimated the mean flux ratio of FeI,XXV and XXVI 
Ka lines as 0.2 : 1.0 : 0.2. Thus the FeXXV Ka fiux is 
about 70% of the total flux of Fe Ka line complex. Taking 
this flux ratio into account, our results are consistent with 
that of Revnivtsev et al. (2006). We, therefore, do not 
discuss origin of FeXXV Ka line in the GRXE further 
than Revnivtsev et al. (2006). 
4-1.2. The combined regions of the GCDX and GRXE 

In the previous section 4.1.1, we found that FeXXV Ka 
profile in the GCDX shows clear excess above the SMD 
model which gave a nice fit to that in the GRXE. As is 
shown in figure 4a, this excess in the GCDX is confined in 
the region of the NSC and NSD. Therefore we fitted the 
whole data of FeXXV Ka line (the GCDX-I-GRXE region) 
simultaneously, fixing equ = egb = 5.85 x lO'^'^ photons 
s~^ Mq and two free parameters of Ensc smd Ensd- The 
best-fit results are given in figures 4c and d, which show 
globally a good fit with the enhanced parameter values 
of Ensc and Ensd- The best-fit Ensc and £nsd are listed 
in table 4. As is found from table 4, we needed 3.8 ±0.3 
and 19 ± 6 times higher emissivities in the NCD and NSC 
region (GCDX) than that in the GD-hGB region (GRXE). 

The fitting in which the mass emissivities of the NSD 
and NSC (ensd and Ensc) are fixed to the same value of 
the GD and GB, but the stellar masses of the NSC and 
NSD (pc and p^) arc free parameters (see equation Al 
and A2) gives that respective stellar masses of the NSC 
and the NSD are 3.8 ± 0.3 and 19 ± 6 times larger than 
those obtained with the NIR observations (Launhardt et 
al. 2002), normalized to the Galactic ridge. 
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In section 4.3, we will discuss the origin of FeXXV Ka 
line in the GCDX based on these results. 

4-2. Systematic error for FeXXV Ka line profile and the 
SMD model 

Since the statistics is limited, present spectral analy- 
sis in the small region (0.°lx0.°2) cannot confine the in- 
terstellar absorption (A^h) within reasonable errors. We 
therefore did not correct the absorption effect for the 
flux of FeXXV Ka. The largest interstellar absorption 
is A^H ^ 6 X lO'^^ cm~^ near the Galactic center (Predehl 
1995). Using the cross section of the photoelectric absorp- 
tion of Morrison & McCammon (1983), we estimated that 
the flux of Fe XXV Ka is reduced only by 9% in the largest 
case of -/Vh = 6 x 10^^ cm~^. Since the largest reduction is 
expected in the GCDX region, the excess discussed in sec- 
tion 4.1.2 becomes even larger by about 10%, if we correct 
the interstellar absorption. A larger error is in the mak- 
ing of the SMD model, which is estimated to be about 
50% (Muno et al. 2006) (see Appendix). Taking these 
errors into account, our conclusion of the large excess of 
Fe XXV Ka line at the GCDX to the GRXE in the best-fit 
SMD model is still valid, and hence no essential change of 
the discussion is required. 

4.3. Origin of FeXXV Ka line in the GCDX 

We found that the mass emissivity of FeXXV Ka line 
in the GCDX region (NSC and NSD) shows excess of 
3.8±0.3 (NSC) and 19±6 (NSD) times larger than that in 
the GRXE region (GD and GB), or the stellar masses in 
the NSC and NSD regions are 3.8 ±0.3 and 19 ±6 times 
larger than those estimated with the NIR observations. 
Integrating over the whole region of the NSC and NSD, 
we obtain that the total stellar mass in the NSC and NSD 
are about 1.4 x 10^ A/q for the former case, which is fully 
consistent with the kinematic mass distribution (see fig- 
ure 14 of Launhardt et al. 2002). However the latter case 
requires total stellar mass within the NSC and NSD to be 
5.3 X 10^ Mq, which is unacceptably large. We hence ig- 
nore the latter possibility, and discuss based on the former 
possibility. 

Origin of Fe XXV Ka line in the GCDX therefore cannot 
be the same as that in the GRXE. Yamauchi et al. (2009) 
found that the line intensity ratio of Fe XXVI Ka to Fe XXV 
Ka in |/*| < 0.°3 is higher than that in |/*| > 8°, which 
also supports that origin of Fe XXV Ka line in the GCDX 
may be different from that in the GRXE. Revnivtsev et 
al. (2009) resolved more than 80% of the Fe Ka flux at 
(/, b) = (0.°1,— 1.°4) into point sources. Although this fact 
supports that origin of FeXXV Ka line in the GRXE is 
due to point sources, that in the GCDX is not yet an- 
swered. As is shown by the arrows in figures 4b and 4d, 
Fe XXV Ka from this region is due to the sum of nearly 
equal contribution of the GD and GB (the GRXE com- 
ponents) and not the GCDX components (the NSC and 
NSD). 

If we assume that origin of Fe XXV Ka line in the GCDX 
is an ensemble of point sources with the same populations 
in the GRXE region, the iron abundances in the GCDX 



should be 3.8(±0.3)-19(±6) times larger than that in the 
GRXE. This is unlikely because the abundance of iron has 
been reported at most two times of solar in the GCDX 
region (e.g. Koyama ct al. 2007b; Nobukawa et al. 2010). 

Another possibility is that the point source population 
in the GCDX is different from that in the GRXE. We 
see the FeXXV Ka intensity excess in the positive /* side 
(east) compared to the negative side (west) at |^*| = 0.°1- 
0.°3 (sec figure 3a). This excess is globally associated to 
the clusters of high mass stars such as the Arches and 
Quintuplet clusters. Muno et al. (2009) also reported that 
the number density of X-ray stars in these regions exceeds 
that in the surroundings. They proposed that the excess 
is due to massive stars associated to the high activity of 
star formation in this region. The most probable mas- 
sive stars are WR or OB stars. These stars often show 
moderately large equivalent width of FeXXV Ka (EWe.r) 
similar to the GCDX (e.g. Skinner et al. 2007; Hyodo 
ct al. 2008). These massive stars might contribute to not 
only the Arches and Quintuplet clusters but also the entire 
GCDX. These stars are relatively bright with the typical 
X-ray luminosity of ^ 10'^''^ erg s~^ (Portegies Zwart et al. 
2002; Wang et al. 2006). The deep Chandra observation 
revealed that point sources with the luminosity larger than 
1032 erg s-i contribute less than 10% of the GCDX (see 
figure 6 of Revnivtsev et al. 2007). We therefore concluded 
that massive stars do not mainly contribute to Fe XXV Ka 
line in the GCDX. 

On the other hand, it may be conceivable that faint X- 
ray point sources with the flux below the current detection 
limits largely contribute to the GCDX. The most probable 
candidate sources for the Fe XXV Ka line are cataclysmic 
variables (CVs), because their spectra are represented by 
high temperature plasma with large EWe.7 and are asso- 
ciated with the 6.4 and 7.0 keV lines, very similar to those 
of the GCDX in the Fe Ka line band. The mean EWqj 
for X-ray bright CVs in the luminosity range of 1030~33 
erg s~^ is about 200 eV (Ezuka & Ishida 1999; Rana et 
al. 2006), while those of the GCDX and GRXE are 450- 
550 eV (Yamauchi et al. 2009). Thus even for the GRXE, 
the faint (unresolved X-ray) CVs should have, at least, 
two times larger EWqj than the bright CVs. In addi- 
tion, the number density of such (X-ray) faint CVs rel- 
ative to near infrared stars should be 3.8(±0.3)-19(±6) 
times larger than that in the GRXE to explain the ob- 
served FeXXV line excess in the GCDX. This possibility 
does not violate the total mass constraint of the NSC and 
NSD (^ 1.4 X 10^ Mq), because the mass contribution of 
CVs to that of the SMD (mass of infrared stars) is only 
a small fraction (see e.g., Pretorius et al. 2007 and refer- 
ences there in). For relatively bright CVs, however, such 
number density excess in the GCDX region has not been 
reported (Revnivtsev et al. 2008). For faint CVs, studies 
for the spatial profile have not been done well. In order 
to settle these potential problems in the point source ori- 
gin, deep observations, with lower detection limit than the 
current value by orders of more than 1-2 are required. 

An alternative possibility is that the majority of the 
GCDX (NSC and NSD) comes from diffuse optically thin 
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Table 4. Best-fit result with the stellar mass distribution 
model. 



Component 



Mass cmissivity of 
the FeXXV Ka line e 



Nuclear stellar cluster 

Nuclear stellar disk 

Galactic disk''' 



111 ± 37 

22 ± 2 
5.85 (fixed) 



* See equation A7. The units arc 10 photons s^ A/q . 
^ Mass emissivities of the Galactic bulge and the Galactic 
disk are fixed to the best-fit of |/*| > l.°5 and |6*| > 0.°5 
regions. 

thermal plasma as already proposed by Koyama et al. 
(2007b). Although the total flux of the FeXXV Ka line 
in the NSC and NSD regions is hardly explained by point 
sources, the spatial distributions are similar to the shapes 
of the NSC and NSD. The diffuse plasma, therefore, 
should be closely related to the ensemble of point sources. 
One possibility is that the plasma is made by multiple su- 
pernova explosions. This may be conceivable because the 
NSC and NSD regions are distinguished from the other 
regions (GB and GD) in respect to supernova activities; 
i.e., very high density of molecular gas and on-going star 
formation. 

5. Summary 

1. We have obtained intensity profiles of the Fel, 
Fe XVV, and Fe XXVI Ka lines in the Galactic center 
region of -3° < / < 2° and -2° < 6 < 1°. 

2. The intensity profile of FeXXV Ka is nicely fitted 
with the SMD model in the GRXE region (|;* | > l.°5 
or |6,| > 0.°6), while that in the GCDX (|/*| < l.°5 
and \b^\ < 0.°6) shows 3.8(±0.3)-19(±6) times excess 
over the best-fit SMD model to the profile of the 
GRXE region. 



Appendix. Stellar Mass 
the Milky Way Galaxy 



Distribution Model of 



This appendix describes the details of the three- 
dimensional stellar mass distribution (SMD) and actual 
two-dimensional fitting model (SMD model) in section 
4. The SMD model was originally compiled by Muno 
et al. (2006) using the results of Launhardt et al. (2002) 
and Kent et al. (1991) based on near infrared (NIR) ob- 
servations with COBE, IRAS and IRT (Boggess et al. 
1992; Clegg 1980; Koch et al. 1982). The major emis- 
sion sources at the NIR wavelength (1-7 fim) are main- 
sequence stars and thus the NIR emission is a good tracer 
of the stellar mass. The stellar mass model consists of the 
components of the nuclear stellar cluster (NSC), the nu- 
clear stellar disk (NSD), Galactic bulge (GB) and Galactic 
disk (GD). In order to describe the three-dimensional stel- 
lar mass distribution of these components, we defined the 
X-Y-Z coordinate, where Sgr A* is located at the origin 
{x,y,z) = (0,0,0), and the Galactic plane is the X-Y plane. 



The sun is assumed to be at (8500, 0, 16) in units of pc. 

The stellar mass density model p* (Mq pc~^) is p, = 
Pnsc + Pnsd + Pgd + Pgb ■ The details of each component 
are described below. 

The NSC component is 

Pc 



Pnsc — 



1 + (r/rc)"-^ ' 



(Al) 



where r = {x + y + z-) '- is the distance from the 
coordinate origin (0, 0, 0), Sgr A*. The parameters 
are r^ = 0.22 pc, pc = 3.3 x 10^ Mg pc"^ (r < 6 pc), 
9.0 X 10^ Mq pc-3 {6<r< 200 pc) and ric = 2 (r < 6 pc), 
3 (6^r< 200 pc). 

The NSD and GD components arc respectively given 



PNSD ~ Pd 



rd 



exp 



2d 



and 



PGD = Pgd cxp 



rgdj 



exp 



^Kd 



(A2) 



(A3) 



where r = (x^ +y^y^^ and z are the distance from Sgr A* 
in the Galactic plane and the height from the Galactic 
plane, respectively. The parameters for the NSD are 
rd = 1 pc, Zd = 45 pc, and n^ = 0.1 (r < 120 pc), 3.5 
(120 ^ r < 220 pc), 10 (220 £ r < 2000 pc). We took 
Pd = 300 Mq pc~'^ for r < 120 pc, and adjusted p^ to 
make the function continuous at r =120 and 220 pc. The 
parameters for the GD are pgd — 5 Mq pc~'^, r^d = 2700 
pc and 2:gd=200 pc. 

The GB component is modeled as a triaxial "generalized 
ellipsoidal" bar (Athanassoula et al. 1990; Freudenreich 
1998). Since the major axis of this "ellipsoid" is tilted 
from any of the X-, Y- and Z-axis, we defined a new co- 
ordinate X'-Y'-Z' by rotating the X-Y-Z coordinate along 
the Z-axis by 15° to the east and along the Y-axis by 1° 
to the north. In the new coordinate X'-Y'-Z', the GB is 
given by; 

(A4) 
(A5) 



PGB =Pgb exp(-rs), 



ri_ 






i/c|i 



l/c 



(A6) 



-3 



a^=1100 pc, 
1.6. 



The parameters are pgb = 8 Mq pc 
aj,=360 pc, az=220 pc. c\\ = 3.2 and cj 

For the three-dimensional mass distribution function 
p{x,y,z) of the components NSC (Al) , NSD (A2), GD 
(A3) and GB (A4~6), we defined the mass emissivity of 
FeXXV Ka line e (in units of photons s~^ -^o^) to be the 
ratio of FeXXV Ka photons emitted from a unit stellar 
(infrared stars) mass per second. The number of FeXXV 
Ka photons emitted per second from a small volume in 
the solid angle fl (a rectangle of 0.°1 x 0.°2) and depth As 
at the distance s from the Sun toward (Z* ,6*) direction 
is given as fls^Asep. The observed flux AF{l^,,b^,) from 
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Fig. 4. (a) The Fexxv Ka line intensity distribution along the Galactic plane as a function of It (distance from Sgr A*). For 
simplicity, only negative Z* region is shown, except for those at |i»| > 8.°0 because the model in these region (Galactic disk) is 
symmetric with respect to the positive and negative U. The green, purple, orange and yellow lines are the best-fit SMD model of 
the nuclear stellar cluster (NSC), nuclear stellar disk (NSD), Galactic disk (GD) and Galactic bulge (GB), respectively. The black 
solid line is the sum of the four components of NSC, NSD, GD and GB. The emissivities of these components are fixed to be the 
same. The data of |Z*| > l.°5 and |6t| > 0.°5 regions arc used simultaneously in the fitting. 

(b) The same as the left panel (a), but those along the Galactic longitude at U = — 0.°114 as a function of 6, (distance from the 
Galactic plane). The arrow shows the position of 6 = — 1.°4 which is almost the same region as Rcvnivtsev ct al. (2009). For 
simplicity, only negative fe» region is shown. 

(c) The same as the upper panel (a) but with the best-fit model in the whole region (GRXE and GCDX) with free parameters of 
the emissivities (or stellar mass densities) of NSC and NSD. 

(d) The same as the upper panel (b) but with the best-fit model same as the panel (c). 



the small volume should be Os^Ase/9/47rs^ = ilepAs/ATT. 
Then the observed flux of Fe XXV Ka is given by the line- 
of-sight integration along s, 



F{lt,bt,) ~ -— i (ensc Pnsc + £nsd Pnsd 
47r Jo 

+eGD Pgd + egb PG^)ds. 



(A7) 



This equation is used for the fitting of observed intensity 
profile of Fe XXV Ka line (figure 4 in the text) . The overall 
uncertainty in building of this model is estimated to be 
about 50% (Muno et al. 2006). 
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